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Abstract 

In this paper, observational constraints on the cosmic equa- 
tion of state of dark energy (p = wp) have been investigated 
using gravitational lensing statistics. A likelihood analysis of 
the lens surveys has been carried out to constrain the cosmolog- 
ical parameters Q m and w. Constraints on fi m and w are ob- 
tained in three different models of galaxy evolution: no evolution 
model (comoving number density of galaxies remain constant), 
Volmerange and Guiderdoni model and fast merging model. The 
last two models consider the number evolution of galaxies in ad- 
dition to the luminosity evolution. The likelihood analysis shows 
that for the no-evolution case w < —0.04 and Q m < 0.90 at la 
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(68% confidence level). Similarly for the Volmerange & Guider- 
doni Model the constraints are w < —0.04 and Q m < 0.91 at 
1(7. In fast merging model the constraint become weaker and 
it allows almost the entire range of parameters. For the case 
of constant A (w = —1), all the models permit Q m = 0.3 with 
68% CL which is consistent with the value of Q m inferred from 
various other cosmological observations. 
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1 Introduction 



Recent observations are in concordance with flat cosmological 
models in which the universe is in an accelerating phase. Anal- 
ysis of the magnitude-redshift data of high-redshift type la su- 
pernovae (SNe la) suggests that the ratio of the matter energy 
to the critical energy, fi m , is ~ 0.3.^^'^ SNe la calibrated 
"standard" candles appear fainter than would be expected if 
the expansion was slowing due to gravity. Recent studies of 
the Cosmic Microwave Background (CMB) data also favour a 
nearly flat universe: the first acoustic peak in the angular power 
spectrum of the CMB is located at I ~ (197 ± 6)M Theoreti- 
cal modeling of structure formation based on cold dark matter 
models (CDM) with Q m = 1 fails to reconcile with the obser- 
vations at the quantitative level. By contrast, a spatially flat A 
CDM universe (non-zero cosmological constant) with Q m ~ 0.3 
can explain observations of galaxy clustering on large scales, in- 
creases the age of the universe (which helps to accommodate 
the age of globular clusters), and makes the total energy density 
equal to critical density as generically predicted by inflationary 
cosmology .IH>H 

Neither observational evidence, nor inflationary considera- 
tions tell us that the cosmological term (dominant, negative 
pressure term) is a constant. There are various phenomeno- 
logical models of dynamical-A (dark energy component) present 
in the literature. These are: 

1. A varies with cosmic time or with the scale factor of Friedmann- 
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Robertson- Walker (FRW) metriciSBH 

2. An X-matter cosmology where this unknown form of energy 
is characterised by an equation of state p x = wp x J^^^^ 

3. Rolling scalar field models (quintessence): The A term is 
considered to be a new, physical, classical field with some 
phenomenological Lagrangian. In this class the most popu- 
lar models are those of an evolving scalar field that couples 
minimally to gravity .SHI EH 

In the present work, we focus our attention on the dark energy 
component characterised by the equation of state, namely the 
ratio w = p/p. In particular, we work with an equation of state 
with — 1 < w < because this range fits current cosmological 
observations best. The cosmological constant is a special case 
of this equation of state corresponding to w = — 1. Constraints 
from large scale structure (LSS) and CMB complemented by 
the SNe la data require 0.6 < Q x < 0.7 and w < —0.6 for a flat 
universe^'^ and w < —0.4 for universes with arbitrary spatial 
curvature. Alcaniz and Lima (1999)0 have derived the limits 
on the cosmic equation of state from age measurements of old 
high redshift galaxies (OHRG). They show that if, as indicated 
from dynamical measurements, the density parameter Q m ~ 0.3, 
then w < —0.27. By combining "cosmic concordance" with 
maximum likelihood estimator, Wang et al. (2000)^ find that 
the best-fit model lies in the range Q m = 0.33 ± 0.05 with an 
effective equation of state w ~ —0.65 ± 0.07. 

In this work we use the statistics of gravitational lensing of 
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quasars to set quantitative limits on the present density of dark 
energy component and matter. 

The main aim of this paper is to constrain the cosmic equa- 
tion of state using gravitational lensing statistics as a tool, tak- 
ing into account the evolution of lenses (galaxies). Many obser- 
vations suggest that the galaxies we see today could have evolved 
from the merging of smaller subsystems. Hence, the inclusion 
of the fact that the number density of lensing galaxies changes 
with time is very important in the lensing statistics. 1301^ We 
consider two different evolutionary models of galaxies : the fast 
merger model of Broadhurst, Ellis & Glazebrook (1992)1 and 
the Volmerange & Guiderdoni model (here after VG model) .0 
Both the models consider the number evolution of galaxies in 
addition to the pure luminosity evolution. 

The paper is organized as follows. In Section 2, we review 
the "dark energy component" model. In Section 3, we discuss 
the two evolutionary models we are considering. In Section 4, 
we write down the lensing probabilities for the three models: 
non-evolutionary model (where the comoving number density of 
lensing galaxies remains constant), fast-merging and VG model. 
We also discuss the likelihood analysis of the lens surveys. Sec- 
tion 5 is devoted to a discussion of the results. 

2 Field Equations and Distance Formula 

For spatially flat, homogeneous, and isotropic cosmologies with 
non-relativistic matter and a separately conserved "dark energy 
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component" with equation of state p = wp, the Einstein field 
equations are: 



(i) 



and 



a 



la a 



(2) 



a 



where the overdot denotes derivatives with respect to time, H 
is the Hubble parameter and a is the scale factor. Subscript o 
refers to present values. Q TO and Q x are present day matter and 
dark energy density parameters. As we are mostly interested 
in effects that occured at redshift z < 5, we neglect the energy 
density of radiation. The condition for the flat universe becomes 
Q x = 1 — Q m . For more details see for example Zhu (2000)0 

The proper distance between a source at redshift z and the 
observer at z = is given by 



where xi z ) i s the comoving distance to the source. 

The angular diameter distance for our models, between red- 
shift zl and Zg reads: 



d(0,z) = a x{z) = cH 1 J q 



z dy 



V^(l + ?/) 3 + (l-^)(l + 2/) 3 ( 1+w ) 

(3) 




(4) 
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3 Evolution Of Galaxies 



Gravitational lens statistics and galaxy evolution are linked with 
each other because the galaxies which are evolving through dif- 
ferent mechanisms are basically acting as lenses. Therefore grav- 
itational lens statistics is certainly going to be affected if the 
lenses evolve. The theory of the formation and evolution of 
galaxies is one of the major unsolved problems of astrophysics. 
According to one view, galaxies evolved through a complex se- 
ries of interactions and have settled in the present day form JH'I^ 
Others believe that galaxies were created in a well defined event 
at a very early time in the history of the universe J^'E^ It re- 
mains unclear which process dominates the formation of ellip- 
tical galaxies. Among the many theories of galaxy formation, 
the idea that galaxies may form by the accumulation of smaller 
star-forming subsystems has recently received much attention. 
Observational evidence also supports this "bottom - up" scheme. 
Deep Hubble Space Telescope (HST) images^ and the size- 
redshift relation of luminous early- type galaxies (E/SO) and mid- 
type spiral galaxies (Sabc) indicate that these objects were as- 
sembled largely before z = 1 and have been evolving passively 
for z < 1. Moreover, HST and ground-based telescopes show 
that the galaxy-merger rate was higher in the past and it roughly 
increases with redshift.B'i'S These arguments suggest that the 
galaxies we see today could have been assembled from the merg- 
ing of smaller systems sometime during z > 1. 

Recent observations^ show that there is a strong deficit of 
galaxies with extreme red colors (seen in different separated 
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fields and at different flux limits) than predicted by models in 
which elliptical galaxies completed their star formation by z ~ 5. 
Therefore, elliptical galaxies must have had significant star for- 
mation at z < 5 through merging and associated star bursts. 
The formation of elliptical galaxies in this way is also consistent 
with the predictions of hierarchical clustering models of galaxy 
formation. 

The other piece of evidence which supports the merging hy- 
pothesis comes from the excess of Faint Blue Galaxies [FBG] 
which have been found in many deep imaging studiesH'^SH! 
Comparison with a model which assumes no luminosity evolu- 
tion in the galaxy population shows that at B = 24, the actual 
galaxy count exceeds the model predictions by a factor of 5. 
Merging of galaxies can also solve the surprisingly steep increase 
in the number density of galaxies. I'll HI It is also found that the 
"FBG problem" cannot be resolved in the conventional scenario 
of formation and evolution of ellipticals. In this scenario, ellip- 
tical galaxies are assumed to have formed in an instantaneous 
burst of star formation at high redshifts with no subsequent star 
formation events. 

Although galaxy mergers are no longer a matter of dispute, 
there is no agreement on the current and past rate of galaxy 
mergers. Several authors have attempted to find the dependence 
of merger rate on the redshift.^^ There are many challenges 
for the mergers theory which may be met in the near future 
with more observations. Recently there have been several semi- 
analytic models, motivated by cosmological theory, which may 
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eventually provide a greater understanding of galaxy formation 
and evolution. IHHH The traditional galaxy number-count mod- 
els, however, are still powerful tools in exploring the formation 
and evolution of galaxies and can be treated as complements to 
the semi- analytic techniques. 

To explain the galaxy number-counts various number-luminosity 
evolution models have been constructed. Along with the no- 
evolution model, we consider two different models of galaxy evo- 
lution which try to explain some of the observational facts listed 
above. These models are: 

1. VG modelB 

2. Fast merging model J 

The general philosophy behind these merger models is to as- 
sume that the current giant galaxies result from the merging of 
a number of smaller "building blocks" . It is found that under 
some general assumptions the theoretical predictions of merging 
models nicely fit the observed galaxy counts. The assumptions 
underlying these models are: 

1. Stability of the Schechter form for the luminosity function. 

2. Conservation of total comoving mass density. 

The evolving luminosity function at any redshift z is given as 
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L*(z) being the characteristic luminosity at the knee, (f)*(z) a 
characteristic density and a is the index of faint end slope. 



No-Evolution Model 

This model assumes that the comoving number density of galax- 
ies is constant and the mass of galaxies does not change with 
cosmic time. Therefore 

(f)*(z) = 0*(O) = constant (6) 

The characteristic luminosity at any redshift remains constant 
and hence the mass of galaxy at any redshift is constant. 

L*(z) = L*(0) = constant (7) 
The "0" refers to present-day values. 

The Schechter form is the most commonly used luminosity func- 
tion for the lens galaxies. Therefore eq.(f5j) becomes 

<5>{L,z)dL = <5>{L,z = 0)dL = 0, (-J exp^-—j — (8) 

where 0*, a and L* are the normalization factor, the index of 
faint-end slope and the characteristic luminosity at the present 
epoch respectively. These values are fixed in order to fit the 
current luminosities and densities of galaxies. This is known as 
the Schechter form of the luminosity function.^H 

VG Model 

In 1990, Volmerange and Guiderdoni,^ proposed a unifying 
model to explain faint galaxy counts as well as observational 
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properties of distant radio galaxies. This model of galaxy evolu- 
tion is based on number evolution in addition to pure luminos- 
ity evolution. According to this model the present day galaxies 
result from the merging of a large number of building blocks 
and the comoving number of these building blocks evolves as 
(1 + z) 1 - 5 . 

It is argued that the present luminosity function is the well 
known Schechter Luminosity Functional given in eq. ()5j) above. 
Then at high z, the galaxies follow a new luminosity function 
where (f)*(z) and L*(z) vary as: 



and 



M?) = 0*(O) (! + £)' 



L«(z) = L*(0)(l + z; 



eq.(|5|) now becomes 
<S>{L,z)dL = (l + z) 2l] 



f (1 + ^ 

.1*1* 



exp 



L 



1 + z)' 



(9) 



(10) 



dL 

(ii) 



or 



z) dL = (1 + zf 71 $(L(1 + z)\ z = 0) dL. (12) 

It is seen that the value 77 = 1.5 gives a fair fit to the data on 
high redshift galaxies. The functional form has the following 
properties: 

(i) Self-similarity as suggested by the Press-Schechter formal- 
ism subject to the constraint that the total comoving mass 
of associated material is conserved.^ 
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(ii) The comoving number density evolves as (fr*(z) = 0*(O) (1+ 
z) 11 and the characteristic luminosity of the galaxy luminos- 
ity function varies as L*(z) = L*(0) (1 + z)~ 1] . Thus the 
galaxies are more numerous and less massive for increasing 
z. 

Fast Merging Model 

The fast merger model of Broadhurst, Ellis & Glazebrook 
(1992) [BEG] was originally motivated by the faint galaxy pop- 
ulation counts. i This model also assumes the comoving number 
density of the lenses to be a function of the look back time St 
or redshift z and hence: 

fa(St) = f(St) 0,(0) (13) 

and 

L*(6t) = \f(6t)]- 1 L< t (0)- (14) 

Since luminosity is related to the velocity dispersion of the dark 
halo of the lensing galaxy through the Faber-Jackson relation 
[L oc t> 7 ),0 this form implies that if we had n galaxies at time 
St each with velocity dispersion v, they would by today have 
merged into one galaxy with a velocity dispersion [f(5t)]~v. The 
strength and the time dependence of merging is described by the 
function f(St): 

f{5t) = exp (Q Hq St) (15) 
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where Hq is the Hubble constant at the present epoch and Q 
represents the merging rate. As suggested by BEG,i we take Q 
= 4. The look back time 5t is related to the redshift z through 



(l + y)- l dy 



Now we can rewrite eq.(O) in terms of St instead of z. 



<$>(L,5t) dL = f(5t) 2 (j), 



L 



W) 



exp 



-L 



W) 



dL 



or 



(16) 



(17) 



<S>{L,5t)dL = f(5t) 2 [$(Lf(6t), z = 0)]dL (18) 



4 Likelihood Analysis of Lens Surveys 

4.1 Basic Equations of Gravitational Lensing Statis- 
tics 



For simplicity we use the Singular Isothermal Model (SIS) for 
the lens mass distribution. The cross-section for lensing events 
for the SIS model is given 

(7 = 16 ^ ( H)4 ( M££ ) » i (19) 
c Dos 

where v is the velocity dispersion of the dark halo of the lensing 
galaxy, Dol is the angular diameter distance to the lens, Dos 
is the angular diameter distance to the source and Dls is the 
angular diameter distance between the lens and the source. The 
mean image separation for the lens at z^ takes a simple form 
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Ae = 8n( V -)^. 



(20) 



c D os 

The differential probability dr of a beam having a lensing 
event in traversing dzL is 



dr = nAz) a — — dzi, 
dz L 



[21) 



where til(z) is comoving number density of the lenses and the 
quantity cdt/dzL is calculated to be 



cdt 



a 



dz L (1 + z L ) v/Q TO (l + z L f + (1 - Q m )(l + z L f^ 



■w 



=, (22) 



where a is scale factor at the present epoch. Substituting for o 
from eq. ([Hp , we get 



c/r = nAz) 



16tt 3 4/ D lD ls , 2 1 

— ^ ( n ) 

a Vos a c 



cHl 



-v 



Y- 



cdt 



dz 



dzi 



L 



(23) 



The No-Evolution Model 

The luminosity function is given by eq. (§) and therefore we 
have 

< n L (z) v 4 >= (1 + z L f vt §{L)dL . (24) 

We assume that v is linked with L through the Faber-Jackson 
relation for elliptical galaxies 



' L 



v 



(25) 
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Hence eg . becomes 



4 

< n L (z) v A >= (1 + z L f vi f §(L)dL ) 7 . (26) 

The optical depth can be written as 

dT = ni + ZL f(^M.) 2 l^L dzLi (27) 

V a L>05 / a dz L 

with 

167T 3 A I 4 \ 

F * = ail ^ [ a + 7 + • (28) 

F* is a dimensionless quantity which governs the probability 
of a beam encountering a lensing object. It is a measure of 
the effectiveness of matter in producing multiple images. S For 
our calculation we use Schechter and lens parameters for E/SO 
galaxies as suggested by J. Loveday et al. (hereafter LP EM 
parameters) :fi 0* = 3.2 ± 0.17/i 3 10 _3 Mpcr 3 , a = 0.2, 7 = 4, 
u* = 205.3 /cm/s and F* = 0.010. 

The differential optical depth of lensing in traversing dzL with 
angular separation between and <j) + d(j), in the presence of 
evolution of galaxies, is given by 



d 2 r 

dzidcj) 



d(f)dz L = F*(l + z L ] 



' DqlDls \ 2 1 cdt 7 /2 
v a Dos ) a dz L T(a + 1 + 



D 



OS 



\D LS 



exp 



l-( 



D 



OS 



\D LS 



]^dz L (29) 
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for ellipticals (lenticulars), where = A6/87r(v*/c) 2 with v* the 
velocity dispersion corresponding to the characteristic luminos- 
ity in ( g5p . 

The Evolutionary Models 

In the case where the comoving number density of the lenses 
changes with time, the equations (|2^), ( |27] ) and ( |29|) read as^ 

4 

< n L (z) v l >= (1 + z L f v 4 , J q z)dL [—j , (30) 

dr = F*(l + ZL f(i\z2^) 2 *(f t )(H) If*^, (31) 

and 



d 2 r 

dZL(i(j) 



d(j)dz L = F* (1 + z L ) 3 (— 
x*(<5f) 2+a 



DolDlsY 1 aft 



2 



OS 



D 



LS 



a dzi 
d(j) 

— dz L , 



(32) 



where ^z(^t) = exp(Q Ho5t) for the fast merging model and 
ty(5t) = (1 + z) h5 for the VG model for evolution of lensing 
galaxies. We notice that when 7 = 4 the differential probability 
is the same for both evolutionary model and no evolutionary 
model. 

Lensing increases the apparent brightness of a quasar causing 
over-representation of multiply imaged quasars in a flux-limited 
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sample. This effect is called the magnification bias. The bias 
factor for a quasar at redshift z with apparent magnitude m is 
given by^^'SSH 

B(m, z) = Ml B{m, z, M , M 2 ), (33) 

where 



w x « fdNnV 1 rM 2 dM dN , „ _ , , N 
B (m> z, M,, M 2 ) = 2 log(M), s). 

(34) 

In the above equation (dNq(m : z) / dm) is the measure of num- 
ber of quasars with magnitudes in the interval (m, m + dm) at 
redshift z. 

We can allow the upper magnification cut off M 2 to be in- 
finite, though in practice we set it to be M 2 = 10 4 . Mq is the 
minimum magnification of a multiply imaged source and for the 
SIS model M = 2. 

We use Kochanek's "best modeF'S for the quasar luminosity 
function: 



^(m, Z) OC (l(T a (") + 10 -Km-m)yl_ (35) 

dm 

where the bright-end slope a and faint-end slope b are constants, 
and the break magnitude m evolves with redshift: 

m + (z — 1) for £ < 1, 
m = < m G for 1 < 2 < 3, (36) 

m D — 0.7(2 — 3) for 2 > 3. 
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Fitting this model to the quasar luminosity function data in 
Hartwich & SchadeS for z > 1, Kochanek finds that "the best 
model" has a = 1.07±0.07, b = 0.27±0.07 and m = 18.92±0.16 
at B magnitude. 

The magnitude corrected probability, for the quasar i with 
apparent magnitude rrii and redshift Zi to get lensed is: 

Pi = r(zi)B(mi,Zi). (37) 

Selection effects are caused by limitations on dynamic range, 
limitations on resolution and presence of confusing sources such 
as stars. Therefore we must include a selection function to cor- 
rect the probabilities. In the SIS model the selection function 
is modeled by the maximum magnitude difference Am(6) that 
can be detected for two images separated by A9. This is equiv- 
alent to a limit on the flux ratio (/ > 1) between two images 
f = io°- 4Am ( e ). The total magnification of images becomes 
Mf = Mo(f + l)/(f — 1). The survey can only detect lenses 
with magnifications larger than Mf. This sets the lower limit 
on the magnification. Therefore the M\ in the bias function 
[7| gets replaced by Mf{9). The corrected lensing probability 
and image separation distribution function for a single source at 
redshift zs are given as^ 



Pi (m,z)-p,J B (m,z,M ,M 2 ) ' W 



and 
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where 



(39) 



Pc (A6) 



1 rzs d 2 T 



dz L 



(40) 



t(z s ) h dz L d{A6) 



Equation ( 1391) defines the configuration probability. It is the 
probability that the lensed quasar i is lensed with the observed 
image separation. 

To get selection function corrected probabilities, we divide 
our sample into two parts - the ground based surveys and the 
HST survey. We use the selection functions as suggested by 
Kochanek.Hl 

In our present calculations we do not consider the extinction 
effects due to the presence of dust in the lensing galaxies. 

The image separation distribution, dN/d(A6), is given bySS 




dN N 
lumber of optical quasars 




Vi(A6) = Ml B(m h M f (0), M 2 ) J" dz L . (42) 




The sum of the lensing probabilities p i for the optical QSOs 
gives the expected number of lensed quasars, Ni ens = E v\- The 
summation is over the given quasar sample. 
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4.2 Testing the models against observations 

We perform maximum-likelihood analysis to determine the con- 
fidence level of w and Q m for the case where comoving number 
density of lenses is constant and for the two evolutionary models, 
namely, the VG model and the fast-merging model. 
The likelihood function is 

N v N L 

C= 11(1- Pi) IlPkPck- (43) 

i=l k=l 

Here Nl is the number of multiple-imaged lensed quasars, Nu 
is the number of unlensed quasars, p k: the probability of quasar 
k to get lensed, is given by eq. (ffif) and p l ckl the configuration 
probability, is given by eq. We considered a total of 867 

(z > 1) high luminosity optical quasars which include 5 lenses. 
These are taken from optical lens surveys such as the HST Snap- 
shot survey ,0 the Crampton survey ,0 the Yee survey ,^Surdej 
survey,^ the NOT Survey^ and the FKS survey^ . The lens 
surveys and quasar catalogs usually use V magnitudes, so we 
transform my to a B-band magnitude using B — V = 0.2 as 
suggested by Bachall et al. (1992)0 

5 Results And Discussion 

The likelihood function as defined in eq. is a function of 

the parameters Q m and w. We allow the parameters Q m and w 
to vary in the range < Q m < 1 and — 1 < w < and find the 
maximum of the likelihood C max . The logarithm of the ratio of 
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the likelihood to its maximum — 2 ln(£/C max ) is asymptotically 
distributed like a x 2 distribution with the degrees of freedom 
equal to the parameters involved. In the figures we present 
likelihood function as a function of two parameters and therefore 
the 68%(lcr) and 95.4% (2a) confidence levels are defined by the 
contours on which C is 31.6% and 4.5% of C max respectively. 
For one parameter fitting, — 2 \n(C/C max ) is distributed like a 
X 2 distribution with one degree of freedom. The 68% and 95.4% 
confidence limits on the parameter are where C is 60.6% and 
14% of C max respectively (Kochanek (1993)^1 and Lampton et 
al. (1976)0) . 

• No-Evolution Model : Fig. 1 shows the contours of constant 
likelihood (95.4% and 68%) in the two parameter space (w, O m ). 
The best fit (C max ) occurs for w = —0.33 and Q TO = 0.0. We 
see that w < —0.04 and Q m < 0.90 at la (68% confidence 
level). For the case of constant A i.e. w = —1.0, C max occurs 
for tt m = 0.44 (fi A = 0.56) and 0.27 < Q m < 0.75 at la. We 
get Q A < 0.83 at 2a for w = -1. Waga & MiceliS have done a 
similar study with Kochaneck's values^ for Schechter and lens 
parameters. For a one parameter fit they get Qa < 0.55 at 2a. 

• VG Model : Fig. 2 summarises the results for the VG Model. 
The maximum of likelihood occurs for w = —0.33 and Q m = 0.0 
with w < —0.04 and Q m < 0.91 at la. For the constant A case, 
C max occurs for O m = 0.44. Also 0.26 < Q m < 0.74 at la and 
tt m > 0.17 at 2a. 

• Fast Merging Model : Fig. 3 shows the contours of constant 
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likelihood (68% and 95.4%) for the Fast Merging Model of galaxy 
evolution. The best fit occurs for w = —0.26 and Q m = 0.0. For 
this model of galaxy evolution the constraints are very weak. For 
a one parameter fit i.e. for the case of constant A, C max occurs 
for n m = 0.52 and 0.29 < tt m < 0.92 at la and Q m > 0.16 at 
2a. 

Table 1 summarises the results for the three models. 

Torres and Waga, i Waga and MicelS have done a similar study 
for decaying A cosmologies (p\ oc a~ m ) without taking evolution 
of galaxies into account. They point out the fact that the con- 
straints obtained from gravitational lens statistics are very weak 
in these decaying A cosmologies. The reason for this is that the 
distance to an object is smaller (for fixed z and Q m ) in decay- 
ing A models as compared to the distance of the same object 
in constant A models for the same value of z and Q m . Hence 
the probability of lensing of that object is reduced in decaying 
A model cosmologies. Therefore the constraints on m from the 
statistics of lensing become weaker. 

In this paper our main aim is to check how the constraints 
on w changes when we include galaxy evolution into the lens- 
ing statistics ( in dynamical A models : X - matter cosmology 
i.e. p x oc a~ 3 ( 1+w ^) . It is a well known fact that the lensing 
probability depends upon the source redshift, luminosity of the 
quasar, lens model, the cosmological model and more impor- 
tantly on the number density of galaxies (lenses). The evolution 
of galaxies changes the number density of galaxies and hence 
the probability of lensing changes. 
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Fig. 4 illustrates how the optical depth (lensing probability) 
decreases with increasing w. This is because as w increases, 
the distance to an object at redshift z decreases which in turn 
decreases the probability of lensing. The inclusion of galaxy 
evolution further decreases the lensing probability in dynamical 
A models ( X - matter cosmology ). As a result the constraints 
on the cosmological parameters (w, Q m ) become weaker. In the 
case of a two parameter fit, we observe that higher values of w 
permit only low values of Q m . 

The probability for a quasar to get lensed decreases marginally 
when we go from no-evolution model to VG model for galaxy 
evolution and hence the constraints on w and Q m obtained in 
the case of the VG model are similar to those obtained in the 
no-evolution model. On the other hand, the probability for a 
quasar to get lensed decreases considerably when we work with 
Fast Merging model for galaxy evolution This gets reflected as 
weaker constraints in the Fast Merging model. Jain et alJ^ con- 
strained f^A (constant A) for flat cosmologies taking into account 
galaxy evolution. They worked with the BFSP quasar sam- 
ple. They constrained it using the information of the number 
of lensed quasars and concluded that evolution permits a larger 
value of Q\. They also found that in their formalism 7 = 4 
masks the difference which the evolution of galaxies makes. We 
would like to point out that the present formalism allows us 
to study the difference which comes in when the evolution of 
galaxies is taken into account even with 7 = 4. 

It is interesting to note that the image separation distribution 
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is different for the three models (Fig. 5). This image separation 
is sensitive to lens, Schechter and cosmological parameters.^ 
Fig. 6 gives the predicted number of lenses Ni ens in the adopted 
quasar sample for different Q m values (w = — 1 case) in the 
three models. The results are similar to those obtained by Jain 
et al. i.e.^the evolutionary models permit larger values of Qa- 
The constraint that in the given optical sample Nl = 5 gives us 
different values of £1a in different models. 

1. No evolution: Q A ~ 0.54 

2. VG model: Q A ~ 0.56 

3. Fast merging model: ft\ ~ 0.60 

In the calculation of number of lensed quasars, Nl in the 
given sample (Fig. 6), no information of image morphology is 
taken into account as explained in Sec. 4.1. Likelihood analy- 
sis, on the other hand, takes into account the observed image 
separations. There is not much difference in the predicted val- 
ues of Qa by these two methods in the case of no evolution and 
the VG model for (w = —1) as can be seen by comparing the 
above mentioned values to those in Table 1. However, for the 
constant A case, likelihood analysis shows that the Fast Merging 
model permits higher values of Q m . Calculating the predicted 
number of lensed quasars and constraining the parameters using 
the observational fact that in the adopted quasar sample there 
are five lenses does not make use of the information about the 
image separations observed for the lensed quasars. The results 
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of likelihood analysis, however, crucially depend on the infor- 
mation about image separation. We feel that it is important to 
incorporate the image separations information while predicting 
(constraining) the cosmological parameters. 

The constraints obtained from lensing statistics are sensitive 
to both lens and Schechter parameters. For the quasar lu- 
minosity function used in this paper, it is observed that for the 
case of one parameter fit ( w = — 1) K96 parameters (Table 2) 
give Nl = 5 ( for this optical sample of 867 quasars) only when 
f2 m is nearly equal to one. For all other values of fi TO , K96 pa- 
rameters give Nl greater than the observed number of lensed 
quasars (there are five observed lensed quasars) in this opti- 
cal quasar sample. The LPEM parameters, on the other hand, 
give numbers around this observed number of lensed quasars 
over a wide range of parameter (Q m ) values. We have therefore 
restricted ourselves to the use of LPEM lens and Schechter pa- 
rameters to study how the galaxy evolution models change the 
constraints on the cosmological paramters. We observe that the 
probability for a quasar to get lensed increases with decrease in 
Q m or/and w. Therefore for small values of Q m and w we can- 
not use the approximation ln(l — p\) ~ —p i (in eq. ^) which 
holds when p\ <C 1. Infact, for consistency, we have not worked 
with this approximation for the entire range of the cosmological 
parameters. 

It is interesting to note that with K96 parameters, the lensing 
probability of a few bright quasars at high redshift becomes 
greater than one for Q m ~ 0.0 and w ~ — 1. For instance, one 
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quasar in the sample with z = 3.02 and tub = 16.1 gives a 
probablity(^) of 1.42 for Q m = 0.0 and w = — 1. This is not an 
artifact of our numerical programming and can be checked easily 
by using either eq. ( pSQ of this paper or eq. (4.15) of K93.^ It 
may be because of the specific values of parameters that we 
have taken. Specifically, the value of F* for K96 parameters 
(F* = 0.026) is quite large as compared to that for the LPEM 
parameters (F* = 0.010). Clearly this larger value of F* pushes 
up the probabilities significantly. In fact, even for those sample 
points where the probability is less than one, there are instances 
where the probability is large enough ( near Q m ~ 0.0 ) so that 
the approximation ln(l — ~ — p i is invalid. 

It could also be that the formalism itself is inconsistent, though 
there is no obvious inconsistency when parameters other than 
K96 are used. If the problem is that of the formalism not being 
consistent, then one would need to investigate exactly what is 
it in the formalism which leads to these seemingly absurd prob- 
ability values. 

The strength of the constraints obtained also depends strongly 
on the lensing data. Extended surveys are required to estab- 
lish n(A6) as a powerful tool for constraining the parameters. 
The upcoming Sloan Digital Sky Survey which is going down to 
1 — a magnitude limit of ~ 23 will definitely increase our under- 
standing of lensing phenomena and cosmological parameters. A 
large number of new gravitational lens systems (~ 18) have been 
discovered by the Cosmic Lens All-Sky Survey (CLASS)S^l . 
CLASS is the largest survey with more than 10,000 radio sources 
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down to a 5 GHz flux density of 30 mJy. But the major disad- 
vantage of a radio lens survey is that there is little information on 
the redshift-dependent number-magnitude relation. This leads 
to serious systematic uncertainties in the derived cosmological 
constraints. 

In this work we have included the effects of galaxy evolu- 
tion to constrain the cosmological parameters (w, Q m ) obtained 
from lensing statistics. We find that while the constraints ob- 
tained in the case of VG model are similar to those obtained 
for the no evolution case, they are much weaker in the case of 
the Fast Merging model. We have only considered those galaxy 
evolutionary models which are based on the assumption that 
the total comoving mass is conserved during mergers. There is 
a possibility that the total comoving mass isn't conserved when 
mergers take place. Work is under progress to see how lensing 
statistics and cosmological parameters are affected when total 
comoving mass density of lenses is not conserved. 

Acknowledgements 

We are thankful to I. Waga for providing us the quasar sam- 
ple. The authors are greatful to C. Kochanek, I. Waga and A. 
Mukherjee for useful discussions during the course of this work. 
One of the author (A. Dev) thanks University Grant Commis- 
sion of India for providing research fellowship. 

References 



27 



[1] Bahcall, J. N., et al. 1992, ApJ, 387, 56 

[2] Baugh, C. M.,Cole, S., k Frenk, C. S. 1996, MNRAS, 282, 
L27 

[3] Baugh, C. M., Cole, S., k Frenk, C. S. 1996, MNRAS, 283, 
1361 

[4] Baugh, C. M., et al. 1998, ApJ, 498, 504 

[5] Bloomfield Torres, L. F., k Waga, I. 1996, MNRAS, 279, 
712 

[6] Broadhurst, T., Ellis, R., k Glazebrook, K. 1992, Nature, 
355, 55 [BEG] 

[7] Burkey, J. M., et al. 1994, ApJ, 429, L13 

[8] Caldwell, R. R., Dave, R., k Steinhardt, P. J. 1998 Phy. 
Rev. Lett, 80, 1582 

[9] Carlberg, R. C, et al. 1994, ApJ, 435, 540 

[10] Chiba, T., Sugiyama, N., k Nakamura, T. 1997, MNRAS, 
289, L5 

[11] Chiba, T., Sugiyama, N., k Nakamura, T. 1998, MNRAS, 
301, 72 

[12] Chiba, M. k Yoshii, Y. 1999, ApJ, 510, 42 

[13] Cole, S., et al. 1994, MNRAS, 271, 781 

[14] Crampton, D., McClure, R. D., k Fletcher, J. M. 1992, 
ApJ, 392, 23 

28 



[15] de Bernardis, P., et al. 2000, Nature, 404, 955 

[16] Driver, S. P., et al. 1995, ApJ, 449, L23 

[17] Efstathiou, G. 1999, preprint astro-ph/9904356 

[18] Eggen, O. J., Lynden-Bell, D., k Sandage, A. R. 1962, ApJ, 
136, 748 

[19] Ellis, R. S., et al. 1996, MNRAS, 280, 235 

[20] Ellis, R. S. 1997, ARAkA, 35, 389 

[21] Faber, S. M., et al. 1976, ApJ, 204, 668 

[22] Fukugita, M., k Turner, E. L. 1991, MNRAS, 253, 99 

[23] Fukugita, M., Futamase, T., Kasai, M., k Turner, E. L. 
1992, ApJ, 393, 3 

[24] Glazebrook, K., et al. 1995, MNRAS, 273, 157 

[25] Guiderdoni, B., k Rocca-Volmerange, B. 1991, A&A, 252, 
435 

[26] Gunn, J. E., k Gott, J. R. 1997, ApJ, 176, 1 

[27] Hartwich, F. D. A., k Schade, D. 1990, ARAkA, 28, 437 

[28] Huterer, D., k Turner, M. S. 1999, Phy. Rev. D, 60, 081301 

[29] Jain, D., Panchapakesan, N., Mahajan, S., k Bhatia, V. B. 
1998, Int.J. Mod. Phys., A13, 4227 

[30] Jain, D., Panchapakesan, N., Mahajan, S., k Bhatia, V. B. 
2000, Mod. Phys. Lett A, 15, 41 

29 



[31] Jain, D. et al., 2001, preprint astro-ph/0105551 

[32] Jaunsen, A. O., Jablonski, M., Petterson, B. R., & Stabell, 
R. 1995, A&A, 300, 323 

[33] Kochanek, C. S. 1993, ApJ, 419, 12 [K93] 

[34] Kochanek, C. S. 1996, ApJ, 466, 638 [K96] 

[35] Kochanek, C. S., Falco, E. E., & Schild, R. 1995, ApJ, 452, 
109 

[36] Lampton, M., Margon, B., k Bowyer, S. 1976, ApJ, 208, 
177 

[37] Lima, J. A. S., & Alcaniz, J. S. 2000, MNRAS, 317, 893 

[38] Lilly, S. J., et al. 1995, ApJ, 455, 108 (1995) 

[39] Loveday, J., Peterson, B. A., Efstathiou, G.,& Maddox, S. 
J. 1992, ApJ, 390, 338 (LPEM) 

[40] Mao, S., & Kochanek, C. S. 1994, MNRAS, 268, 569 

[41] Maoz, D., et al. 1993, ApJ, 409, 28 (Snapshot) 

[42] Ostriker, P.J., & Steinhardt, P. J. 1995, Nature, 377, 600 

[43] Partridge, R. B., k Peebles, P. J. E. 1967, ApJ, 147, 868 

[44] Peebles, P. J. E. 1984, ApJ, 284, 439 

[45] Perlmutter, S., et al. 1999, ApJ, 517, 565 

[46] Perlmutter, S., et al. 1999 Phy. Rev. Lett, 83, 670 (1999) 

30 



[47] Press, W. H., & Schechter, P. 1974, ApJ, 187, 487 

[48] Ratra, B., k Peebles, P. J. E. 1988 Phy. Rev. D, 37, 3406 

[49] Riess, A. G., et al. 1998, A J., 116, 1009 

[50] Rix, H. W., Maoz, D., Turner, E., & Fukugita, M. 1994, 
ApJ, 435, 49 

[51] Rocca-Volmerange, B., & Guiderdoni, B. 1990, MNRAS, 
247, 166 

[52] Rusin, D & Tegmark, M. Preprint No. astro - 
ph/0008329 (2000) 

[53] Sahni, V., & Starobinsky, A. 2000, Int. J. Mod. Phys., D9, 
373 

[54] Schechter, P. 1976, ApJ, 203, 297 

[55] Schwezier, F. 1996, A J, 111, 109 

[56] Silveira, V., k Waga, I. 1994, Phys. Rev. D, 50, 4890 

[57] Silveira, V., & Waga, I. 1997, Phys. Rev. D, 56, 4625 

[58] Surdej, J., et al. 1993, A J, 105, 2064 

[59] Takahashi, R. & Chiba, T. Preprint No. astro - 
ph/0106176 (2001) 

[60] Toomre, A. 1997, in The Evolution of Galaxies and Stellar 
Populations eds: B. M. Tinsley & R. B. Larson ( Yale Univ. 
Observatory), p-401 

31 



[61] Turner, E. L.,Ostriker, J. P., & Gott III, J. R. 1984, ApJ, 
284, 1 

[62] Turner, E. L. 1990, ApJ, 365, L43 

[63] Turner, M. S., & White, M. 1997, Phy. Rev. D, 56, 4439 

[64] Waga, I., & Miceli, A. P. M. R. 1999, Phy. Rev. D, 59, 
103507 

[65] Wang, et al. 2000, ApJ, 530, 17 

[66] Yee, H. K. C, Filippenko, A. V., k Tang, D. 1993, A J, 105, 
7 

[67] Yee, H. K. C, & Ellingson, E. 1995, ApJ, 445, 37 
[68] Zepf, S. E., & Koo, D. C. 1998, ApJ, 337, 34 
[69] Zepf, S. E. 1997, Nature, 390, 377 
[70] Zhu, Z. H. 2000, Int. J. Mod. Phys., D9, 591 



32 



Model 


Fit 


Best Fit 


68% CL 


95.4% CL 






w, tt m 






No-Evolution 


A 


-0.33, 0.0 


w < -0.04 










tt m < 0.90 






B 


0.44 


0.27 < tt m < 0.75 


0.17 < tt m 


VG 


A 


-0.33, 0.0 


w < -0.04 










tt m < 0.91 






B 


0.44 


0.26 < tt m < 0.74 


0.17 < fl m 


Fast Merging 


A 


-0.26, 0.0 








B 


0.52 


0.29 < Vt m < 0.92 


0.16 < tt m 



Table 1 : Here A refers to the two parameter fit and B refers to one parameter 
fit i.e. when w — — 1. 



Survey 


a 


(f)*(h 3 Mpc- 3 ) 


v*(km/s) 


p* 


LPEM 


+0.2 


3.2 x 10- :i 


205.3 


0.010 


K96 


-1.0 


6.1 x 10~ 3 


225.0 


0.026 



Table 2: Lens and Schecter Paramters for E/SO galaxies 
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No Evolution 
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Figure 1: Contours of constant likelihood (68% and 95.4%) for No-Evolution 
Model for lensing galaxies 
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VG Model 

10 1 1 1 1 1 1 — 
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Figure 2: Contours of constant likelihood (68% and 95.4%) for Volmerange 
& Guiderdoni Model of galaxy evolution. 
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Fast - Merging 
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Figure 3: Contours of constant likelihood (68% and 95.4%) for Fast Merging 
Model of galaxy evolution. 
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Figure 4: Probability of lensing as a function of w 
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AG ( arcsec) 



Figure 5: Predicted image-separation distribution dN/d(A9) compared with 
the observed image-separation in the optical quasar sample (histogram). 
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Figure 6: Pridicted number of lensed quasars with A8 < 4" in the adopted 
optical quasars sample. 
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